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1. INTRODUCTION 

DC motors are commonly used in industry applications [1], [2]. A significant difficulty for DC-motor 
control design is the unknown time-varying nature of its parameters [3], [4]. An important framework for motor 
control design is the state adaptive backstepping (SAB) technique presented in [5], as can be noticed from [6], [7]. 
Nussbaum gain techniques are usually incorporated in the SAB control framework in order to handle the effect 
of unknown time-varying parameters and to improve the robustness of the system. Robust SAB control schemes 
incorporate a compensation term in the control law, and some modification in the update law, for instance the c 
modification, see [8]. Nevertheless, upper or lower bounds of the plant model coefficients have to be known to 
guarantee asymptotic convergence of the tracking error to a residual set of user-defined size. In [9] a non-adaptive 
state backstepping control scheme is developed. The resulting time derivative of the Lyapunov function involves 
an unknown, time varying but bounded term, whose upper bound is unknown, such that the backstepping states 
remain bounded and converge to residual set whose size depends on plant parameters. 


Nussbaum SAB control schemes are based on the controllers presented in [10], [11], as can be noticed 
from [12], [13], [14]. In turn, the controllers in [10], [11] are based on the Universal Stabilizer that was originally 
introduced in [15] and discussed in [16]. The main drawback of the Nussbaum gain technique is that the result- 
ing upper bound of the transient behavior of the tracking error depends on integral terms that involve Nussbaum 
functions [12], [17], [14], [10]. Some controllers that use this technique usually present some of the following 
drawbacks: i) some upper or lower bounds of plant model parameters have to be known in order to guarantee the 
convergence of the tracking error to a residual set of user-defined size [14], [10], but the control designs in [12], 
[13] indicate that a proper design would relax this drawback, and ii) the control or update laws involve signum type 
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signals, as can be noticed from [17]. 


In [18], a linear induction motor is considered. The friction force and unknown time-varying model 
parameters lead to a lumped bounded uncertain term whose upper bound is unknown. The goal is to control the 
mover position. The drawback is that the identification error is assumed constant in the definition of the Lyapunov 
function. In [6], a synchronous motor driven through AC/DC rectifiers and DC/AC inverters is considered. It is 
assumed that the motor parameters experience unknown time-varying but bounded behavior. The goal is to control 
the motor speed, the rectifier output voltage and the d component of the stator current; the tracking error converges 
to a residual set whose size depends on unknown motor parameters and user-defined controller parameters. There- 
fore, if upper bounds of unknown motor parameters are known and controller parameters are properly chosen, the 
size of the residual set can be user-defined also. 

Other works address the problem of designing a controller for different motors [19, 20, 21]; despite the 
fact that the controlled systems operate as it is expected, the main disadvantages of these works are: the size of 
the output error cannot be determined, and no analyses of the system behavior inlcuding measurement noise are 
presented. In [22] two coupled controllers are designed: a Linear Quadratic Gaussian (LQG) and a MRAS-based 
Learning Feed-Forward Controller (LFFC); eventhough the simulation results demonstrate the potential benefits of 
the proposed controlled, the LQG algorithm may fail to ensure closed-loop stability when variations in the uncer- 
tainties are large enough. In [23] a plant model in controllable form with unknown varying but bounded parameters 
is considered, being the upper bounds of such parameters unknown. A SAB control scheme is developed, and over- 
comes the main drawback of the Nussbaum gain method, as a result, the transient error is upper bounded by an 
unknown constant that does not depend on integral terms. Nevertheless, the control scheme is only valid for plant 
models in “companion form". The method is based on the Lyapunov-like function technique appearing in [24], 
[25]. 


In addition to the undesired unknown time-varying nature of plant parameters, other important issue is 
the measurement noise. Tracking performance can be degraded, even if the controller is robust against modeling 
uncertainty and disturbances (cf. [26], [27]). Some of the main techniques to tackle the effect of measurement 
noise are: high gain observers, interval observers, filter theory and the technique developed in [9]. High gain 
Observers are useful to estimate system states and output derivatives (see [28], [29]). In [28], a nonlinear plant 
model in state-space form and a plant model in controllable form, are considered, respectively. Both plant models 
involve known constant coefficients. The real output is defined as the first state, and is measured, whereas the other 
states are not. The output measurement is expressed as the sum of the real output plus a bounded measurement 
noise parameter. The observer depends on the difference between the noisy output measurement and the output 
estimate. The stability analysis indicates that the state estimation error converges to a residual set whose upper 
bound depends on the magnitude of both the measurement noise and the observer parameters. Such upper bound 
has a global minimum for some value of the observer gain (see [28] and [29]). The main drawbacks of the design 
are: the size of the residual set is unknown, so that the upper limit of the steady state of the state estimation error is 
unknown, and second, the coefficients of the plant model are required to be known. Interval observers provide an 
upper and a lower bound for each unmeasured state variable (see [30], [31]). The main disadvantage of the interval 
Observers is that several upper and lower bounds of the plant model parameters are required to be known. 


In [9], it is assumed that the measurements of the plant states are corrupted by noise and are described 
by a measurement model, which consists of a polynomial with respect to the real state vector, with degree one 
and unknown but bounded time-varying coefficients, being the time derivatives of such coefficients unknown and 
bounded. If each measurement model is differentiated with respect to time, the time derivative of the state mea- 
surement is a linear polynomial of degree one with respect to the time derivative of the real state vector. The states 
resulting from the backstepping state transformation are defined in terms of the noisy measurements instead of 
the real states. To compute the time derivative of each quadratic function of the backstepping procedure, each 
measurement model is differentiated. In [32], the state adaptive backstepping (SAB) of [5] is used as the basic 
framework for controlling a DC permanent magnet motor whose voltage is supplied by a buck power converter, 
but the effect of noisy measurement is not taken into account. 


In the present paper, significant modifications are incorporated to the control scheme presented in [32], 
in order to handle the effect of measurement noise: i) a measurement model is used to define the relationship 
between the motor current (ia) and the motor speed (W), and their corresponding noisy measurements (7,),,, and 
Wmjm, ii) the states of the backstepping state transformation are defined in terms of the noisy measurements, and 
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iii) the differentiation of each quadratic-like function with respect to time involves the differentiation of each mea- 
surement model. The controller is implemented in a digital platform to carry out the real experiments. The main 
contribution of this paper respect to close ones are: a) the relationship between the real and the measured state is 
taken into account in the control design procedure: in the definition and differentiation of the backstepping states 
as well as in the definition and differentiation of the quadratic functions, b) upper or lower bounds of the noise 
model parameters or their combination are not required to be known by the controller, and c) the convergence of 
the tracking error to a residual set of user-defined size is proven in presence of noisy measurments. 


The rest of the paper is organized as follows. In section 2., the plant model and the control goal are 
described. In section 3. the controller is designed. In section 4., the bounded nature of the closed loop signals and 
the convergence of the tracking error are proven. In section 5., numerical and experimental results are presented. 
Finally, discussion and conclusions are presented in Section 6.. 


2. PLANT MODEL AND CONTROL GOAL 
The DC motor is represented by the following plant model (see [33]): 


W. Bw k Thrict Tr 
m — J. m + JÉ ta Cs J. ) 

mM tg ic (1) 
ig = — TL.a A T. m+ L u 


The state variables are: the armature current 7, and the motor speed Wm. u is the control input and it corresponds 
to a voltage value. The system output is y = Wm. The model parameters are: the voltage constant ke [V/rad/s], 
the armature inductance La [H], the armature resistance Ra [Q], the viscous friction coefficient B [N.m/rad/s], the 
inertia moment Jeg [kg.m?], the motor torque constant k; [N.m/A], the friction torque Tfric [N.m], and the load 
torque 77, [N.m]. 

The following assumptions are made for the model (1): Ai) the parameters Tz and Jeq are time-varying, 
unknown and upper bounded by unknown constants, and Jeq is positive and lower bounded by an unknown positive 
constant, Aii) parameters B, kt, Ra, La, ke are unknown, positive and constant, and Aiii) Wm and ia are measured, 
but their measurements Wm|m, tajm are noisy and satisfy: 


Wa = Wim + a7 and talm = Agta + ag (2) 


where the parameters ag, a7, as, and ag are unknown, time-varying and upper bounded by unknown positive 
constants, their time derivatives are unknown, time-varying and bounded by unknown positive constants and the 
parameters ag and ag are positive and lower bounded by unknown positive constants. The above expressions are 
based on [9]. For a simpler control design, the plant model (1) is rewritten as: 


Ly = —a121 + a313 — a3 (3) 
d = —4T] — A5X% + bu (4) 
zı = Wm, T2 la, U = Ve, Y= Wm (5) 
Ym = Wim = agt; + az (6) 
Ləm = iam = Agl2 + ag (7) 
where 
B ky 
ay =, a2 — —— (8) 
Jeq eq 
(Tae + Tr) ke Ra 1 
a3 = 0. — 7, a= —, ag — V—, b= — (9) 
Jeq La La La 


so that a1, a2, à3, G4, a5, and b are positive, a4, a», and a3 are time-varying, a4, as and b are constant, and ym and 
X2m are the noisy measurements of Wm and ta, respectively. Assumption Ai implies that the parameters a1, a2, 
a3 are unknown and time -varying, but they are upper bounded by unknown constants. Assumption Aii implies 
that a2 > 0, and the parameters a4, a5, and b are unknown and constant. Assumption Aiii implies that the states 
zı = y and x2 are unknown but their measurements Ym and x2 are known, the parameters ag, a7, ag, à9, a6, à, 
dg, and àg are unknown and time-varying but bounded, and the parameters ag and ag are positive. In summary, the 
system (1) satisfies the following properties: Pi) a1, a2, a3, a4, a5, b, ag, and ag are positive, Pii) the parameters 
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Q1, Q2, Q3, Q6, A7, Ag, Ag, Q6, Q7, Gg, and åg are unknown, time-varying and upper bounded by unknown constants, 
Piii) the parameters a», ag, ag, and 1/ag are positive and lower bounded by unknown positive constants, Piv) the 
parameters a4, a5, and b are unknown, positive and constant, and Pv) the values of y = x, = Wm and x2 = i, are 
unknown, whereas their measurements ym and £2m are known. 


Taking into account the model, the control goal can be defined as follows. Let the following reference 
model: 


Yd = Am1Yd AmoVd H GmoV mre f (10) 


where W,,,..; is the user-defined reference value, and amı and amo are user-defined positive constants. Hence, 
the desired output y4 is provided by (10) subject to i) amo and amı are user-defined but positive and constant, and 
ii) Wmref is user-defined and bounded but it may be time-varying. Therefore, equation (10) is a stable reference 
model. 

The tracking error is defined as: 


e(t) = ym(t) — yalt) = Winjm — Ya (11) 
Qe = {e: |e] € Coe} (12) 


where ym is defined in (6), ya is the desired output and is provided by (10), Qe is a residual set whose size is 
defined by C, which is an user-defined positive constant. The goal of the control design is to formulate a control 
law and an update law for the plant model (1), subject to assumptions Ai to Aiii, such that: CGi) the tracking error 
e asymptotically converges to the residual set Qe, CGii) the control and update laws do not involve discontinuous 
signals, CGiii) the control law and the updated parameters are bounded, and CGiv) all the closed loop signals are 
bounded. 


3. CONTROL DESIGN 

In this section, a controller for the plant (1) is developed, it takes into account the assumptions Ai to Aiii, 
and the goals CGi to CGiv stated previously. The procedure is similar to that in [32], but there are several differ- 
ences due to the presence of measurement noise. Therefore, the procedure omits the steps that are quite similar to 
those in [32]. The state adaptive backstepping (SAB) presented in [5] is used as control framework, but important 
modifications are incorporated in order to tackle the effect of unknown time-varying plant model coefficients and 
measurement noise. 


The controller design procedure is organized in the following steps: i) define the first new state zı and 
differentiate it with respect to time, ii) define a quadratic function V,; that depends on 24, and differentiate it with 
respect to time, iii) express the terms that involve time-varying coefficients, as functions of upper constant bounds, 
and parameterize such bounds as function of parameter and regression vectors, iv) express the parameter vector 
in terms of updating error vector and update parameter vector, and define the second new state z2, v) differentiate 
z2 with respect to time, define a quadratic function V, that depends on zı and z2, and differentiate it with respect 
to time, vi) express the terms that involve time varying coefficients as function of upper constant bounds, and 
parameterize in terms of parameter and regression vectors, vii) express the parameter vector in terms of updating 
error vector and updated parameter vector, and formulate the control laws, and viii) formulate the Lyapunov-like 
function, differentiate it with respect to time and formulate the update laws. 


Step 0. In this step, the model (3)-(4) is expressed in terms of ym, £2m, using the noise models (6)-(7). 
Differentiating (6) and using (3), yields: 


Ym = (àg = à601)21 + 0206412 — 0306 + a7 (13) 


solving (6) and (7) for xı and x2 and substituting into the above expression we obtain the basic expression for Ym: 


(4g — asaı)  , a2ao a2069 
Um — m 1 Tam G306 
a6 ag ag 
az 
+å7 — (dg — agai). (14) 
6 


Differentiating (7) with respect to time and incorporating (4), yields: 


Loam = —a48X1 + (dg = à58)X2 + agbu + ag (15) 
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solving (6) and (7) for xı and x2 and incorporating in the above equation we obtain the basic expression for tom 
which is 


: a4ag ag — a5ag a7 
Lm = Um + Lam d agbu + agag — 
ag ag a6 
1 ag n 
-F(—àg + à5ag) — + dg (16) 
ag 


Remark 3..1 In order to consider the noisy Ym and X2m instead of the real but unknown values xı and zs, in the 
remaining procedure equations (14) and (16) are used instead of equations (3) and (4). 


Step 1. In this step, the first state variable is defined and differentiated with respect to time. The state 
variable zı is defined as the tracking error: 


Zi — € — Um — Yd (17) 
where yq is provided by (10). Differentiating (17) with respect to time and using (14), yields: 
21 = Ym — Ya (18) 
A= — + E Lam + (a1a6 — D 
TET ase + år — ġa (19) 


Step 2. In this step, a quadratic function that depends on 2; is defined and differentiated with respect to 
time. Such quadratic form is defined as: 


Va = 0/25 (20) 
Differentiating (20) with respect to time, using (19) and adding and substracting c; 27 yields: 
. : aga —ajag +a EN 
Va-unác-adcea--cra + 1 [tes + (aras — àg) — 
ag a6 a6 
G20689 


(21) 


d36 + a7 + C121 — Ya 
ag 

The term —c422 has been added to obtain asymptotic convergence of the tracking error later. The unknown and 

time varying behavior of the bounded parameters a1, a2, a3, Q6, G6, A7, 47, Gg, and ag is a significant obstacle 

for the controller design. For this reason, the terms that involve such coefficients will be expressed as function of 

upper constant bounds. 


Step 3. Recall that a1, a2, a3, a6, G7, Gg, Q9, dg, and åy are time-varying, unknown and bounded. In this 
step, the terms that involve such time-varying parameters are expressed as function of upper and lower constant 
bounds, and such bounds are parameterized in terms of parameter and regression vectors. The term that involves 
the brackets in (21) yields: 


—aj1a6 +46 — , . a7  d2069 ; : 
Zu Um + (a1a6 — àg) azas + G7 + C121 — Ya 
ag ag ag 
< iolymllzi] + gailzi] + lerz1 — yall Zr (22) 


where j11o, H11 are unknown positive constants such that 


Lla1aetáe | < n 
— ae | $510 

a 

° a agaga (23) 
(aias — a6) a e 4 7 2396 + G7] S iai 


As mentioned in [32], Young's inequality must be applied to (22), such that the |z;| term leads to z7, to allow a 
proper definition of z2. Arranging (22), and applying Young's inequality ([34]), yields: 


—a,a6 + 6 . a7 020609 »^ : 
zı | — — — Ym + (a1a6 — 46) 436 + àz + C121 — Ya 
a6 a6 ag 
io m lez — 9allal 
< 6s — Win ||21| 6a ea] + — —— —— — (24) 
Ca Ca Ca 
2 72 2 2 2 2,2 
Co |, Bir 2, Ca , Hio 2 2, Ca , (Cia — 9a) zi 
+ apt F zi 4 + 25 
at EET EE 2d (53) 
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where c, is a positive constant that should be chosen to fulfill certain conditions that will be defined later. Substi- 
tuting into (21), yields: 


2 
. C Q2 
puit p Ei qui 
2 a 2c 


2 
7 2c 


1 10 
2 t2 
Ca Ca 


252 
C121 — Ud) £ 

Ca 
The terms that involve x2), 27, 227, (c121 — Ya)? should be grouped in a new known state variable z2, according 
to the procedure in [5]. Nevertheless, the unknown time-varying behavior of parameter azas /ag poses a significant 
obstacle. To remedy such situation, a positive constant lower bound of azas /ag will be used. Property Pi mentions 
that a2, ag, and ag are positive, whereas property Piii implies that a2, ag, and 1/ag are lower bounded by unknown 
positive constants. Therefore, 


26 


0 < in» € (27) 


where jij12 is an unknown positive constant lower bound. Incorporating the constant /ij12 into (26) and arranging, 
yields: 


i2, p, 
Q206 2 
Va <- Tum Ca + RI BO DI 
zl cz? A — ds T2m 2d Tin au, Paz dy. 2d Tm a, Pam 
1 1 2 28 
23 fi aacra — da) (28) 
cesi. came id 29 
= -C12, + 2! zi ia Zam + P Voge tert (29) 


where q is a known regression vector and 0; is an unknown constant parameter vector given by: 
2 S 2T 
Yi = [1, Um: (ez Ya) ] (30) 


= E T 
is B iio 1 
1 | Se 5 
Hü2 Hi2 H2 


(31) 


Step 4. In this step, the unknown parameter vector 0, is expressed in terms of an updating error vector 
and an updated parameter vector, and then a new state variable z2 is defined. The parameter vector 0, can be 
rewritten as: 


6, = 6, — 0 (32) 
where 
TE -2 -2 1 T 
(opi Eu Sum. ue (33) 
Hi2 Mu2  Hn2 


where Ó, is an updated parameter vector provided by an updating law that will be defined later, and 6, is an 
updating error. Substituting (32) into (29), yields: 


: 3c aga 
2 2%6 
Va < C121 + z + VAI 
2 ag 


~ 1 _ 2:1. 
Lam + 94 01 2 Buzz? — 916i 2 Bz? 
a a 


As can be noticed from [23] and [24] the updated parameter Ó, is non-negative, so that là; | = 6, for 6; (to) > 0. 
The aecompisameut of this property will be shown later. In view of this fact and incorporating the inequality (27) 
in the term yy 16,0. 5c, ? 111227 and arranging, yields: 


; 3c? 206 ; 1 a2a6 
Va < 74+ — 4 T 2.010 34 
z1 € —02j 2 zl us Lam + Q1 32 aa zi — P11 xod (34) 
3c? aja ~ JL _ 
= ez? + 8 t a ome — 9] fie? (35) 
2 ag 2c2 
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where 


; 1 
z 
Z2 = Zam + 916177571 (36) 
2c2 
Step 5. In this step, the state variable 2» is differentiated with respect to time, a quadratic function V, is 
defined as function of zı and z2, and such function is differentiated with respect to time. Differentiating (36) with 
respect to time, yields: 


a 1 Ans ^ 1 
e eee ST T paT : 
22 = Lam + 94101 22^ + yı 01 2d zı «0 23 zl (37) 
where 
G1 = (0, 2ymým, 2(c121 — Sa) (ei — dja)] " (38) 


Substituting (18) into (37) and arranging, yields: 
Z9 = Lam + Y1bYm + Pic (39) 


where 


1b (2 (vo Pius ta(aa- jfi) " 


BEP! 
2c2 


+97 61) (40) 


; : SAN 1 
Pre = —2(c121 — Ya)(C1Ya + Ya) 91/3} a2 


à l a 1 
T T. : 
Tv, 0 2d z—v6 LL (41) 


a 


(1p and q. are known scalar functions, us and is are the second and third entries of the vector 61, and Ya, Ya» 
and 94 are provided by (10). Substituting (14) and (16) into (39) and arranging, yields: 


. _ aag _ (dg — asas) , (a6 — asa) , 4246 
Z2 = Um d Ləm T P1bYm 1 Pibt2m 
ag ag a6 ag 
42060. ] . a7 
ap (-55m — azas + à; — (às — cot) ) 1b 
ag a6 
QT n ag P 
a5 («es + (—ås + agag) — + às) t agbu + pıc (42) 
a6 ag 
The quadratic form that depends on z, and z» is defined as: 
V, = (1/2)(2t + 23) (43) 


Differentiating with respect to time, incorporating (35) and (42), and adding and subtracting —c»22, yields: 


Vz = 2121 + 2922 = Vz + 2020 


2 2,932 
X —c121 — €2% + zC 


2 a 
426 a4ag ag — A5ag ag — aga 
+22 on Um — Tm + — (Q1ibUm 
a a6 
a2a6 20609 : 2 a7 
+ p1pL2m + | —————- — agas + àz — (ås — aga1)— | Yio 
ag ag a6 


+ (aaas SZ + Cri + asas) "2 + tn) + caza + asbu + | 
408 8 T Q508 r ag 222 8 Pic 
a6 ag 


NET 
2257 gi nazi (44) 
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The control input u is defined as follows: 
U = Ua + Up (45) 
where ua is an user-defined constant. The constant ua was incorporated in order to avoid abrupt control input 


behavior at the beginning of the closed loop operation. Therefore, it should be chosen as the open loop value of 
the control input. The control u, is established by the controller design. Substituting (45) into (44), yields: 


3 
2 
Vz € -e1 2] — C223 + gla 
a2a6 a4ag ag — A5ag ag — 601 
+22 (21 + Q1bt2m) — Ym + ———— fam + ———— — (ibm 
ag ag ag a6 
020609 LA . a7 
T G306 T a7 (ås a604) (?1b 
ag ag 


+ | asa S ud dada) + ag | + agbu 4 + coz 
8 a6 8 TO Me a a silo 222 


~ 1 
agbuyz2 — p1 ĝa zgn (46) 


Step 6. Recall that a1, a2, a3, Q6, a7, Ag, A9, å6, å7, åg are unknown and time-varying. In this step, 
the terms that involve such parameters are expressed in terms of upper bounds, and parameterized in terms of 
parameter and regression vectors. The term that involves the squared brackets can be rewritten as: 


az2a6 a4ag ag — a5ag ag — ded 
22 | - (z1 + q1bt2m) a Me Ləm + — —— — ibm 
8 6 8 6 
G20609 . . a7 
ee — a3dg + à; — (dg — a601) — | Pip 
ag ag 


a f a : 
T (mast + (—ås + asag) 2 ES às) + agbua + qc + exa] 
a6 ag 


< |29| [Bas ym] + Baa|x2m | + Bas iem] + Baclzi + q1bt25 
+ft17|Y10| + jas + jao|ua] + eie + c222]] (47) 


where ji13, H14, H15, B16, [417, H18, H19 are unknown positive constant upper bounds that satisfy: 


aaa = ag—a5a = 
a | S Bs, |A | S Aa 
àg—aedi| < q azas | < m 
S Hi5, S Hi6 
E a ia fs a Py (48) 
— 2208 — azae + a7 — (dg — aga1) 27 S lar 
asas + (—àg + asas) + às < fis, |asb| < pio 


From equations (46), (47), and as can be inferred from [32] there are two significant obstacles. First, the |z2| term 
may lead to discontinuous signals in the definition of uy. This can be remedied by using Young's inequality. Sec- 
ond, the unknown varying behavior of agb makes it difficult for u to eliminate the effect of the terms that involve 
unknown time varying parameters. This can be remedied by incorporating a positive constant lower bound of agb. 


From properties Piii and Piv (see page 2182) it follows that ag and b are positive, b is constant and ag is 
lower bounded by an unknown positive constant. Therefore, 


0 < jii29 < aab (49) 


where jio is an unknown positive constant lower bound. Rewritting (47) in terms of a parameter vector and a 
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regression vector, and incorporating (49), yields: 


d206 aag åg — A5dg , a6 — a601 
22 (21 + Q1b22m) Ym d Lom + £1bUm 
dg ag ag ag 


(20609 : I a7 
aN — agag + à; — (dg — aga1)— | 10 
ag a6 


a a 
+ (ua + (—ås + asas) => + às) + agbua + Yie + exa] 
6 8 


< vy Taiz0|z2|8' 02 (50) 
where 
E = [lym]; |E2ml; o1]; [21 4-015220], len, 1, ltal, [Pre + c2z2l]" (51) 
TN _ z _ _ _ 7 
05 = i [H13. faa; Bass Pie, fiz, fis, Ao, 1", (52) 
20 


— is a regression vector whose entries are known, and 05 is a parameter vector, whose entries are positive, constant 
and unknown, and j12o is an unknown positive constant lower bound. The constant ./jij29 has been incorporated 
in order to handle the effect of the unknown time-varying parameter agb appearing in the term agbup in (46). 


Step 7. Since the parameter vector 05 is unknown, in this step it is expressed in terms of an updated 
parameter vector and an updating error vector; after of this the control law is formulated. The parameter 05 can be 
rewritten as 


05 = 02 — 0» (53) 


where 65 is an updated parameter vector provided by an update law which is defined in the step 8, and 05 is an 
updating error vector given by 


1 


0, = 0, — — —][jns, fia, Bis. Bae; fav, fas, Big, 1" (54) 
L120 
Substituting (53) into (50) yields: 
aga asa åg — asa ag — aga 
al 2^9 (21 + qist2m) Ë Um -F ——— 9 — 42m $e ibm 
ag a6 ag ag 


(20609 ; I a7 
ti- — azas + à; — (às — agai) > 1b 
6 


ag 
a : a . 
(«s T4 (—àg + aaa) o + às) + agbtta + qc + cx] 
a6 ag 
< Vii20|22|6 "62 — Viiz0|2216 | 6» (55) 


incorporating the inequality (49) and applying Young's inequality (cf. [34] pp. 123) to the term y Hii20|Z2|G ' 65, 
yields: 


426 a4ag ag — a5ag ag — 601 
z2 (21 + Q15X2m) Ym + Lom + P1bYm 
ag a6 ag ag 
020609 : 3 a7 
+ (- — à3a6 + à; — (àg — aaa) ) (1b 
ag a6 
a f a . 
+ («is T (—àg + a5ag) 2 às) + agbua + qc + = 
ag ag 
c 1 2/ -TÂ 2 = =Tp 
Zu 2,425022 (9 03)* — V üizo|z2]g 02 (56) 


where c, is a positive constant that satisfies some some conditions that will be defined in the Step 8. Substituting 
(56) into (46) and arranging yields: 


; 3 2 1 ^ 
VE c122 C322 | x t E + agbz» (w= za) 
"RE! » 
pı ða 2 F202? — v/lizo|z2| | 02 (57) 
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with the aim to cancel the effect of the term agbz5(1/2)c5 ? z; (i  05)?, the expression for up is chosen as: 


1 -TÂ 2 
Ub = -zal 03) (58) 
In view of (45), the control law for u is: 
1 DATA 
U = Ua — za AG 03) (59) 
Ce 
substituting (59) into (57), yields: 
; , 32,6 TO l- 2 T 516 
V, < —2min(c;, c3] V; 4 aca tt > 1 0, 22 Du202; — VHiz0|22|6 02 


The above expression implies that the time derivative of the Lyapunov function would contain the term (3/ 2)c2 + 
(1/2)c2, so that the required negativeness properties would be altered. Therefore, the quadratic function V, is 
considered, which is a truncated function of V, and vanishes when V, is lower or equal than the constant Cy,;. 
The quadratic function V, is defined as: 


U (1/2)(v Vz — V 6mm if V: > Covz 
V, = ; (60) 
0 otherwise 
Cow = (1/2)CÈ (61) 
where V, is defined in (43). The function defined by (60) and (61) has the following properties: V, > 0, V; < 


3Co- + 3V. and V, and OV, /9V. are locally Lipschitz continuous. Differentiating (60) with respect to time, 
yields: 


uu ATA S (62) 


dt OV, 
where 
U 1VVz-VCbvz ; 
Ve _ J p s EV Os T 
OV, 0 otherwise 
Combining (62) with (60) yields: 
dV. s. OV, 3 c2\ OV, 
V. 2 c z 
dt e a + (s oo 2 
,0V. OV: 
—1 ium ad li2021 == aV, v/lizo| za] ' 0» av. (64) 


Step 8. In this step, the Lyapunov-like function is formulated and differentiated with respect to time, and 
the update laws are formulated. The Lyapunov-like function is defined as: 


Vi) V. Vo (65) 
z(t) = [n (t), a(t), 65, 62] EM (66) 
Vo = (1/2)jt208] Dy 61 + (1/2)v/jir208 T3 62 (67) 


where 6, and A, are defined in (33) and (54) respectively, and V. is defined in (60). The vector Z(t) contains the 
closed loop states z1(t), z2(t), 01, 93 . For the sake of simplicity, V(z(t)) is represented as V. Differentiating 
(65) and (67) with respect to time, yields: 


V - V. Ys (68) 
Vo = jia001 T1 104 + jiio04 T3 165 (69) 
Incorporating (64) and (69) into (68), yields: 


BV. . 3 OV, _ x 1 20V ep 
V« —2min(c, c5) V; oV. G Ê+ ce) + mol ( Piza By, FI å) 


- OV, 
evasit ( - ales zti) (70) 
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To cancel the effect of the terms that involve the updating errors 0, and 65, the update laws are formulated as: 


1 Oir ss 
= HV ioe 2^19y. av,” 04 (to) > 0 (71) 


6, = 
OV. 
ó =f B 72 
2 = l'3lzsle OV. (72) 
where I; and Ig are 3 x 3 and 8 x 8 diagonal matrices, respectively, whose elements are user-defined positive 


constants. The involved state dependent parameters are: pı (30), OV. /OVz (63), @ (51), and z (36). From (30) 
and (63) it follows that 


OV, 80V, 
yı = lal, By v [os . (73) 


(73) together with (71) leads to 6, = LN > 0. Substituting (71) and (72) into (70), yields: 


. OV, 3:5 2 
V< oV. (2min{er,co}¥ - PI :- 3) (74) 


Although the control law (59) and the update laws (71) and (72) have been formulated, the values of the constants 
Ca and ce have not been defined. Expression (74) can be rewritten as: 


OV: 1 3c +c 
V< V,+ Vz - E 75 
~min{e1, c2} avr = ( + oy (75) 

The constants ca and ce are defined as user-defined positive constants that satisfy: 
3c? + c? € 2min{c1, c2}Coyz (76) 


A simple choice that satisfies the above requirement is: 


Ca = c, = (1/2)Cbe y min{ c1, c2} (77) 
From (63) it follows that: 


8V./0V, = 0 if Ve < Cy; 


0V,/0V, > 0 if V, > Cy (78) 


From (78), (76) and (75) it follows that: 
V. 
V < —min{c,, olay Ve (79) 


Figure 1 schematizes the plant model and the proposed controller. 


Remark 3..2 The developed controller involves the control law (59) and the update laws (71) and (72). The state 
dependent signals and the corresponding equations are: zi (17), z2 (36), pı (30), Yio (40), vic (41), v D), 
V, (43), V; (60), and OV, /OV:z (63). The signals ya, Ya, and Yq are provided by equation (10). The user-defined 
positive constants are: Cy, co, the diagonal Glens of 1 and I», Ca, Co Coes Amor and amı; Wmreg is the desired 
reference value. In addition, Cruz = (1/2)C? be» Ca and ce must satisfy (76), and ua is the value of the control input 
u in open loop operation previous to the closed loop operation and it is constant and user-defined. 


Remark 3..3 The user defined constants may be set as follows. i) amo and am; are positive and they may be fixed 
to values such that the desired output ya provided by (10) tracks the refence value W,ep with the transient response 
required by the user. ii) The diagonal elements of I1 and I matrixes are positive and constant (li; > 0); they 


affect the dynamics of 6, and bo, therefore, they must be chosen by trial and error simulations so as to obtain 
a suitable transient behavior of 61, 62 and tracking error e, and a suitable value of the control input u. Small 
values for elements Tij; lead to a slow convergence of tracking error. iii) c, and ce are positive and are set such 
that de? + c2 < 2min(ci, c21 y; according to equation (76); a simple values that satisfy such requirement is 
Ca = Co = V;Cye / min(ei, c2}. iv) Coe is the size of the residual set Qe to which the tracking error converges; 
hence, it should be chosen as the maximun absolut value of the allowed tracking error in steady state. v) cı and co 
are positive constants used in the backstepping procedure. They may be chosen by trial and error so as to acieve 
suitable transient behavior of the tracking error. 
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Wmref Reference] Ya, Yds Ya 
model 


Ud; Ya, Ya 


61,42 


Control law 


0 = sa model m 


Figure 1. Schematic diagram of the plant and the controller. 


4. BOUNDEDNESS AND CONVERGENCE PROPERTIES 
This section presents the bounded nature of the closed loop signals and the convergence property of the 
output. 


Theorem 4..1 Boundedness of the closed loop signals. Consider: i) the plant model given by (1), which is subject 
(43), V. (60), and 0V,/8V, (63), and iii) the signals Ya, Ya, Ya that are provided by equation (10), where amo 
and amı are user-defined positive constants and Wmregș is an user-defined value; and iv) the user-defined positive 
constants c, C2, the diagonal elements of Tı and V2, Ca, Ce Cbe, amo, m1, where cq and ce must satisfy (76). In 
addition, ua is an user-defined constant, and Cy; = (1/2)C7.. If the controller defined in (59), (71), and (72) is 
applied, then the signals z1, z2, Ó,, Ô>, and u remain bounded, and the tracking error e remains upper bounded by 
an unknown constant: 


lel € 2 {Veo + (1/2) i20 A Gy (to) + Tus (to) + UV (to) 


&Q/2) ps [Api 9) 4 eng 9] ) + Che (80) 
where 
pis { (1/2) (V Vzs — VObwz) if Vo 2 Ciz d 
0 otherwise 
Vio = (1/2) Ca (to)? + 22(to)*) (82) 
Voo = (1/2)fir20(O1(to) — 61) "Py (Ôa (to) — 61) 
+(1/2)vV irzo (62(to) — 02) T5 | (62(to) — 02) (83) 
Proof 4..1 From (65) and (79) it follows that 
V «0, V(a(t)) € V(z(t,)) (84) 


where &(to) is the vector Z(t) at t = to, obtained by combining (66) with t = to, whereas V(Z(to)) is the 
Lyapunov-like function at initial time to, obtained by evaluating (65) at t = ty: 


V (x(t5)) z Ve Bi Veo (85) 
s M RE E 
Vio = (1/2) (21 (to)? + z2(to)*) (87) 
Voo = (1/2)ji120(61(to) — 61) "Py (Ôi (to) — 61) 

+(1/2) V irz (82 (to) — 02) T3 (82(to) — 62) (88) 


From (65), (85) and (84) it follows that 
V, < Vo F Voo and Vo < Tos T Veo (89) 
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From the above expression it follows that V, € Læ. To establish the bounded nature of zı and z», a relationship 
between V., 21, and zz has to be established and incorporated in expression (89). Solving (60) for Vz yields: 


= (VI + Us) iV. (90) 


Vz € Cy; otherwise 


Combining (43), (89), (90), yields: 


lal < V2 (V27; + VCs) (91) 
lal € V2(V2W Vio + Voo + V/ Ctuz) (92) 
|za| € V2(V2/ Vio + Voo + V/ Ctuz) (93) 


According to previous results zy € Lo; and z3 € Læ. To prove that 6, and 65 are bounded, the first task is to 
define a notation for each one of the diagonal elements of Tı and V2: 


a = Tip, ip %1/2 = * Tw. 2p 1/3 = zin 3] 

YLS = Ton, 1]: ‘Y2/2 T = 5 Doo, 2] 12/3 = = S Tot, 3] (94) 
2/4 = * Tala, 4) 72/5 = " 5 Dat,s] ^2/e = D'ate,g] 

^óarr Ê Vor,;z] "roja = T2j8,8] 

From the above notations, and equations (67) and (89) it follows that 


= (1/2) 20 Eros if] top ip] + 1139113 | 


(1/2) Riso [nA 8p + + Izba] < Veo + Voo (95) 


The above expression yields: 
(1/2) i0 A811 < Vzo + Veo 


(1/2) Bao; 58s < Vo + Voo 


(96) 
(1/2 ) 11207219511 2[1] as Vo + Voo 


(1/2) v i058 2 < Vis + Voo 


Therefore, ĝi € Læ and 02 € Loo. This result together with definitions (33) and (54) completes the proof for 


ĝi € La and 05 € Lex 


The bounded nature of the control input u is proven at the following. Equations (17) and (30) together 
with zı € Læ, yield xı € Lo, Ym € Loo and pı € Loo. Equations (7) and (36), together with pı € Læ, 
01 € Loo, z1 € Læ, and zg € Ls, yield Lam € Lo and x? € Lo. Equations (43), (63), and (71) jointly with 
zı € Læ, 22 € Læ, and yı € Ls, yield Vz € Ls, 9V,/8V, € Leo, and A, € Loo. Equations (40) and (41) 
together with Ym € Leo, 6, € Loo, 21 € Loo, Yı € Loo, and 6, € Loo, yield pip € Loo and pic € Loo. Equation 
(51) together with Ym € Loo, Tam € Les, Yib € Loo, Z1 € Læ, Pie € Loo, and z2 € Lo, yields p € Læ 
Equation (59) together with z2 € Loo, € Loo, and 05 € Loo, yields u € Loo. This completes the proof for 
u E€ Lo 


Hereafter, the upper bound for the tracking error e is established. From (88) and (94) it follows that the 
quadratic form Vọo can be rewritten as: 


Voo = (1/2) fizo [rA bnlte) + o (to) + 13a (to)| 


+(1/2) Rss [n3 8 (to) + + n (66)| < Veo + Voo (97) 
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From (17), (61), and (91), it follows that |e| < v2( 2V; + (1/V2)Cte), incorporating (89) and arranging, 
yields |e| < (2v Vzo + Voo + Cre), now incorporating (97), we obtain 


lel < 2{Veo + (1/2)fu2o [4/182 (to) + 12a lto) + 75708 lto) 
aa = 1/2 
t (1/2)V/pi20 ha 8 (to) + + Y29218] (to)] ) + Coe (98) 


Notice that expression (80) indicates that the upper bound for e(t) can be rendered small by choosing 
large values of 71/1, Y1/2, Y1/3» Y2/1> ^^^» Y2/8 and a low value of Cre. However, hardware considerations must 
be taken into account. 


Theorem 4..2 Convergence of the tracking error. Consider: i) the plant model given by (1), which is subject to 


V, (60) and OV, /OV. (63), iii) the signals ya, Ya, a provided by equation (10), where amo, amı are user-defined 
positive constants, and Wmref is a user-defined reference value, and iv) the user-defined positive constants c4, 
C2, the diagonal elements of lı and V», Ca, Ce, Coe, Gmo, ami, Where cq and ce must satisfy (76). In addition, 
Cus (1/2). and ua is an user-defined constant. If the controller defined in (59), (71), and (72) is applied, 
then the tracking error e asymptotically converges to a residual set Qe, whose size Cre is user-defined: 


Qe = {e : |e] € Coe} (99) 
Proof 4..2 From (63), (79) it follows that: 
V < =c fa < —af, < OYE > to (100) 
where 
ET (1/2) (/ V: = V Couz)V Vz if V. 2 Covz 
fa = { 0 otherwise (101) 


0 otherwise 


= 2; > 
* 24 GIT = VO) VS Cry T 
Arranging and integrating (100) it is obtained: 


V(z(t)) + af fgdr < V(Z(to)) (103) 


where X(t) is defined in (66), and z(t,) is Z(t) evaluated at t = to. From (103) it follows that f; € Ly. In order 
to apply Barbalat’s Lemma it is necessary to prove that f, € Lo. and Í € Læ. Since V, € La it follows from 
(102) that f, € Loo. Differentiating (102) with respect to time yields: 


oo OF 
fo = ea V (104) 
fo | ee if Ve > Ow. dm 
OV, 0 otherwise 


Notice that O f, /OV. is continuous with respect to V.. Since V, € Loo then Of, /0V. € Loo. Because zı € Ls, 
22 € Læ, t1 € Loo, 2 € Loo, and u € L it follows from (19) and (37) that 21 € Loo and 22 € Los. Thus, from 
(44) it follows that V, € Loo. Because 0f,/0V, € Loo and V, € Loo then it follows from (104) that Ía eu. 
Because fg € Loo and Js € La the Barbalat's Lemma (cf. [35] pp. 76) indicates that f, asymptotically converges 
to zero. From (102), it follows that V; converges to Qyz, where Quz = {Vz : Vz € Cy. Furthermore, from (43) 
it follows that zı asymptotically converges to 1,4, where Qa = {21 : |z| € V2Chuz}. Since Cy,; = (1/2)C7, 
and zı = e, then e asymptotically converges to Qe, where 


Qe = fe: je] € Coe} (106) 
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Remark 4..1 In Theorem 4.2 the asymptotic convergence of the tracking error e to the residual set Qe is proven. 
The size of the residual set is given by Cye which is user-defined, positive and constant; even more, Cre does not 
depend on i) plant model parameters. ii) bounds of plant model parameters. iii) measurement model parameters. 
iv) bounds of measurement model parameters. These features of Cre are achieved by the control design procedure 
and te main tasks to this end are: i) to incorporate unknown positive constant bounds to handle unknown parame- 
ters from the plant and the measurement models. ii) to incorporate these bounds in the parametrizations. iii) to 
incorporate a truncation in the quadratic form Vz. 


5. NUMERICAL AND EXPERIMENTAL RESULTS 

The permanent magnet DC motor (PMDC) has the following characteristics: a rated power of 250 Watts, 
a power voltage of 42 VDC, a rated current of 6 Amps, and a maximum speed of 4000 rpm. The motor speed Wm 
is recorded by an encoder of 1000 pulses per turn. The armature current 7, is measured by a series resistance. The 
backstepping control technique is implemented in the control and development card dSPACE DS1104. This card 
is programmed from Matlab@®)/Simulink platform and it has a graphical display interface called ControlDesk. The 
sampling rate for all variables (Wm and ta) is set to 4 kHz. The state variable 7, is 12-bits resolution; the controlled 
variable Wm is sensed by an encoder which has 28-bits resolution. At each sampling time (250 jus) the controller 
uses the measured Wmm and tajm to calculate the voltage to feed the motor. The control law and the update laws 
defined in section 3., and the parameters Cy, c4, c», 1, l2, Ua and Wmreg are entered to the control block by the 
user. The motor parameters and the user defined controller parameters are stated in Table 1, including the size of 
the residual set Qe, i.e. Cbe, which has been fixed to 5. 


The measured values W,,,,, and tajm are used by the backstepping controller to compute z1, 22, 91, 
Ø, Q1b. Pic» 0V,/8V,, V., and V,, and then to compute ĝi, bo, and u, according to remark 3..2. To test the 
performance of the backstepping controller, simulations and experiments are considered. In the simulation case, 
the motor model (1) is used to generate the values of W,njm and tajm. In this case, the motor parameters presented 
in Table 1 are introduced by the user, in order to allow the numerical simulation of model (1). In the experimental 
case, the values of Wm|m and tajm are measured by sensors. Recall that the values of motor parameters Ra, La, 


B, Jeq» Kt, ke, T fric» and Ty, are not involved in the control and update laws. 


[^] 
e 
e 


280 


(model output) (rad/s) 
N 
3 
(measured output) (rad/s) 
N 
3 


240; 240 
220; 220 
E 200 200 
E J E : j 
0,5 1 15 18 = 05 1 15 18 
Time (s) Time (s) 

Figure 2. Simulated results for step change in Figure 3. Experimental results for step change in 
Wet from 200 rad/s to 300 rad/s att = 1s Wmref from 200 rad/s to 300 rad/s at t = 1 s 


The initial tuning of adaptive parameters can be made off line, on line or using the simulation model. In this case 
off line tuning was made. The reference signal Wmrep is changed from 200 rad/s to 300 rad/s at t = 1 s in closed 
loop, in both the experimental setup and simulations. The transient of W,,),, is shown in Figures 2 and 3; the 
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Motor parameter 

Ra: Armature resistance 

La: Armature inductance 

B: Viscous friction coefficient 
Jeq: Moment of inertia 

k: Motor torque constant 

ke: Voltage Constant 

Tfric: Friction torque 

T: Load torque 

Wer: Reference speed 

F'c: Switching frequency 

F's: Sampling frequency 

1T' p: Unit time delay 

Wm: Motor speed 

ia: Motor current 
Zero-order-hold: Fixed step discrete time 


Value 

2.7289 Q 

1.17 mH 

0.000138 N.m/rad/s 
0.000115 kg.m? 
0.0663 N.m/A 
0.0663 V/rad/s 
0.0284 N.m 
variable N.m 

200 rad/s 

4 kHz 

4 kHz 

250 us 

28 bits quantization 
12 bits quantization 
250 us 


User-defined controller parameter 

Tı: updating gain 

(being 75 the unity matrix of dimension 3) 
T2: updating gain 

(being Is the unity matrix of dimension 8) 
Che 

c, and co: user-defined gains 

Ca and ce: user-defined gains 


Value 
0.0006 J3 


0.0006 Ig 
5 


10 
7.9 


2195 


Ua: initial control action, defined by the user 20 V 
Qm 1: parameter of the reference model (10) 40 
amo: parameter of the reference model (10) — 400 


transient of the tracking error is shown in Figures 4 and 5; and the control input u = ve is shown in Figures 6 
and 7. The steady state error for the simulation and experimental cases is lower than 2.5 96 when Wmre f = 200 
rad/s, and lower than 1.7 % when Wmre f = 300 rad/s. Notice that the user-defined size Cy, = 5 stated in Table 
1 is achieved by both simulation and experiments as the tracking error converges to a residual set of lower size, as 
can be seen in figure 1. During the transient state the error gets below —5 or above 5 because the controller was 
designed to operate in continuous time, whereas the implemented controller works in discrete time; however, in 
steady state the error converges to a pre-defined residual set. 


6. CONCLUSIONS 

The time-varying nature of the parameters, noise and plant models leads to lumped parameter terms, and 
each of those lumped terms has to be treated as an unique parameter. The main tasks to handle the effect of such 
terms are: i) to incorporate unknown positive constant upper bounds, ii) to incorporate unknown positive constant 
lower bounds, iii) to incorporate the lower bounds in the parametrizations, and iv) to incorporate a truncation in 
the quadratic form that depends on the backstepping states. In addition, the application of the Young’s inequal- 
ity allows to avoid the presence of discontinuous signals in the control and update laws. The truncation of the 
quadratic form allows to obtain an adequate time derivative of the Lyapunov-like function, guaranteeing adequate 
boundedness and convergence properties. 


The stability analysis indicates that the transient value of the tracking error is upper bounded by a positive 
constant that depends on initial values of both the system states and updated parameters. This guarantees that the 
transient tracking error can be rendered small by properly choosing: 1) the gains of the update law, and ii) the 
initial value of the desired output. Such bounded nature of the tracking error is an important benefit with respect 
to closely related controllers that use the Nussbaum gain method. 
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Figure 5. Experimental tracking error for step 


Figure 4. Simulation tracking error for step change change in Wmref from 200 rad/s to 300 rad/s at 


in Winref from 200 rad/s to 300 rad/s at t = 1 s t=1s 
40 — 40 
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Figure 6. Simulated control input u — v, for step 
change in W,,,,; from 200 rad/s to 300 rad/s at 
t=1s 


Figure 7. Experimental control input u = v, for 
step change in W,,,.« y from 200 rad/s to 300 rad/s 
att=1s 


The convergence analysis indicates that the residual set for the convergence of the real but unknown error 
depends on the magnitude of the noise parameters. In addition, if the noise parameters are such that noise is absent, 
then the residual set of the real but unknown error is the same as the residual set of the current tracking error. 
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